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Oh' Abstract 

I give an introduction to the theory of meson-antimeson mixing, aiming at students who plan 
to work at a flavour physics experiment or intend to do associated theoretical studies. I derive 
, the formulae for the time evolution of a neutral meson system and show how the mass and width 

Q>^ ' differences among the neutral meson eigenstates and the CP phase in mixing are calculated in 

, the Standard Model. Special emphasis is laid on CP violation, which is covered in detail for 

00 ■ K—K mixing, Bd—Bd mixing and Bg—Bs mixing. I explain the constraints on the apex (p, 77) 

of the unitarity triangle implied by eK, AMb^, AMbj/ AMb^ and various mixing-induced CP 
■"sj" ' asymmetries such as acp{Bd — > J/ipKshoit){t)- The impact of a future measurement of CP 

, violation in flavour-specific Bd decays is also shown. 
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The neutral K, D, and Bg mesons are the only hadrons which mix with their antiparticles. These 
meson states are flavour eigenstates and the corresponding antimesons K, D, B^ and Bs have oppo- 
site flavour quantum numbers: 

K ~ 'sd, D ~ cu, Bd ~ hd, Bg ~ bs, 

K ^ sd, L> ~ cu, Bd ~ bd, Bs ~ bs, (1) 

Here for example ~ bs" means that the Bg meson has the same flavour quantum numbers as the 
quark pair (b, s), i.e. the beauty and strangeness quantum numbers are S = 1 and S = —1, respec- 
tively. The meson states in Eq. (O are also eigenstates of the strong and electromagnetic interactions. 
As long as we neglect the weak interaction, they are also mass eigenstates, with the same mass for 
meson and antimeson. In the Standard Model (SM) all interaction vertices conserve flavour, except 
for the couplings of W bosons to fermionsjl The piece of the SM Lagrangian which describes the W 
couplings to quarks reads 

^^ = ^ E [yjkUjLl''dkLW^ + V*,,dkLl''u,LW-]. (2) 

Here is the weak coupling constant and V is the 3x3 unitary Cabibbo-Kobayashi-Maskawa 
( CKM) matrix: 

(yud Vub \ 

Vcd Vcs Vdo . (3) 
Vu Yts Vth I 

In Eq. ^ I have further used the notations (^1,^2,^3) = [d, s, b) and (ui, U2, M3) = (u, c, t). The 
W boson only couples to the left-handed components of the quark fields as indicated by the subscript 
"L" in Eq. At fourth order in the weak coupling we can change the flavour quantum numbers 
by two units and obtain transitions between mesons and antimesons. The corresponding Feynman 
diagrams are shown in Fig.[T] These |AF| = 2 diagrams, where F denotes the appropriate flavour 

^Strictly speaking, this statement assumes that the so-called unitary gauge for the weak gauge bosons is adopted. The 
unphysical charged pseudo-Goldstone bosons, which appear in other gauges, also have flavour-changing vertices. Changing 
the gauge shuffles terms between the pseudo-Goldstone bosons and the longitudinal components of the gauge bosons. 



1 . 1 Mesons, quarks and box diagrams 



3 




Figure 1: Box diagrams for K—K , D—D , B^—Bd and Bg—Bs mixing. The zigzag lines represent 
W bosons. For each process there is also a second box diagram, obtained by a 90° rotation. 



quantum number F = S, C or B, represent the lowest non-vanishing contribution to the transition 
matrix element S12 defined by 

X4.{4). ^v {M ipM)\S\M jp^)) 

with the S-matrix S and the generic notation M = K,D,Bd or Bg. (The notation E12 refers to the 
quantum-mechanical two-state system with |1) = |M) and |2) = |M).) I comply with the standard 
relativistic normalisation of the meson states, {M{p')\M{p)) = 2E {2'k)^6^^\p' — p). The meson 
mass Mm = yBP^-fP' in the denominator in Eq. (jj]) is introduced for later convenience. In terms 
of the Hamiltonian (density) H^^f^ (x) = — which encodes all interactions of the SM, the 
S-matrix is given by the usual time-ordered exponential 

In order to link Eqs. Q and to the diagrams of Fig.[T]we must consider the contribution from Cw in 
Eq. ^ to —H^J^^ and expand the time-ordered exponential in Eq. ^ to order g^. The determination 
of this term amounts to the calculation of the two contributing box diagrams with the usual Feynman 
rules of the weak interaction. To this point we have only used standard text-book quantum field 
theory, noting an important omission: No effect of the strong interaction has been taken into account 
by now. Most importantly, we do not know yet how to take care of quark confinement, which forces 
the external quarks in the diagrams of Fig. [T] to form mesons. As an important feature. Quantum 
Chromodynamics (QCD) behaves very differently at short and long distances: At short distances 
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(probed by large energies) the QCD coupling constant gs is small and we can apply perturbation 
theory [1], just as we did with the weak interaction. That is, effects of short-distance QCD can be 
included by adding gluons to the diagrams in Fig.[T] At large distances, corresponding to low energies, 
QCD is non-perturbative and one must resort to different methods, such as lattice gauge theory or QCD 
sum rules. Long-distance QCD is also referred to as hadronic physics, because its degrees of freedom 
are hadrons rather than quarks and gluons. In many cases the associated theoretical uncertainties 
are the main obstacle in the relation between measured quantities and the fundamental parameters 
of nature encoded in the Lagrangian C. Theorists pursue a two-fold strategy to deal with hadronic 
uncertainties: On one hand they try to refine non-perturbative methods such as lattice gauge theory. 
On the other hand they try to identify quantities in which hadronic uncertainties are small or even 
absent or look for ways to eliminate hadronic uncertainties through clever combinations of different 
observables. We will encounter both strategies in our discussion of meson-antimeson mixing. Weak 
processes of hadrons involve several largely-separated energy scalesH For example, mB — B mixing 
we encounter nit > My/ » mfc » Aqcd> where Aqcd ~ 0.4 GeV is the fundamental scale of the 
strong interaction governing e.g. the size of binding energies. In order to correctly calculate S12 we 
must separate the different scales from each other and apply different computational methods to large 
and small energy scales. However, without detailed understanding of the strong interaction we can 
roughly assess the relative importance of the contributions from the different internal quark flavours 
in Fig. [T] In the case of Bii — B^ mixing and Bs — Bg mixing one finds that the box diagram with 
internal top quarks vastly dominates over the diagrams with light quarks, because the result of the 
diagram grows with the internal quark mass. For K — K mixing and D — D mixing no such estimate 
is possible, because the contribution with the heaviest quark is suppressed by small CKM elements. 

Owing to S12 / 0, M and M mix and are no more mass eigenstates. The latter are obtained by 
diagonalising the 2x2 matrix Sj^ , where 



with |l,]3i) = \M{pi)) and \2,p2) = \M{p2)) generalises Eq. (01). We list two important aspects of 
meson-antimeson mixing: 

i) The two mass eigenstates are linear combinations of M and M. The degeneracy is lifted and 
we can denote the two mass eigenstates by Mh and Ml, where "H" and "L" stand for "heavy" 
and "light", respectively. Mh and Ml not only differ in their masses, but also in their lifetimes. 

ii) If we produce a meson M at some time t = 0, the corresponding state will evolve into a 
superposition of M and M at later times t > 0. One observes meson-antimeson oscillations. 

We will calculate the differences among the masses and decay widths in the second and third lectures. 
Studies of neutral Kaons mainly exploit property i), while the mixings of the other three neutral meson 
systems are investigated through property ii). The reason for the Kaon's special role here is the vast 
lifetime difference between Kh and Kl- The former state, usually denoted as K\ong, lives roughly 
500 times longer than Kl = i^^short> so that one can easily produce a Kiong beam. For D, B^ and Bg 
mesons the width differences are much smaller than the average decay width of the two eigenstates 
and this method is not feasible. The identification of the meson (discriminating between M and M) 

^^I use natural (or Planck) units with ft = c = 1, so that masses and momenta have units of GeV. 



i{2^f5^'\p[-p,)^. 



(6) 



2Mm 
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needed to track the meson-antimeson oscillations is called flavour tagging. To observe the oscillations 
the mesons must move sufficiently fast in the detector. Modem B factories, which produce {B^, Bd) 
pairs via the T(4S') resonance, have therefore asymmetric beam energies, so that the centre-of-mass 
frame (coinciding with the rest frame of the T(4S')) moves with respect to the laboratory frame. At 
hadron colliders studies of meson-antimeson oscillations profit from the large boost of the produced 
mesons. Tevatron and LHC are especially powerful for Bg physics, because the Bg — Bs oscillations 
are very rapid. 



1.2 A bit of history 



Meson-antimeson mixings belong to the class of flavour-changing neutral current (FCNC) processes, 
which involve different flavours with the same electric charge. Since in the SM such processes are 
forbidden at tree-level, they are sensitive to new heavy particles appearing as virtual particles in loop 
diagrams. Historically, the first new particle predicted from the consideration of FCNCs was the charm 
quark, which was needed to eliminate large tree-level FCNC couplings in conflict with experiment [2]. 
Subsequently, the rough size of the charm quark mass nic was predicted from the size of the mass 
difference AMk = Mh — Ml in the neutral Kaon system [3]. A great success story of flavour 
physics has been the exploration of the discrete symmetries charge conjugation (C), parity (P) and 
time reversal (T). Charged Kaon decays had revealed in 1956 that P and C are not conserved by the 
weak interaction, while physicists kept their faith in a good CP symmetry. If CP were conserved, we 
could assign CP quantum numbers to i^iong and -fiTshort- The latter meson was observed to decay into 
a two-pion state, and each pion is CP-odd and contributes a factor of — 1 to the total CP quantum 
number (which is multiplicative). A further contribution to the CP quantum number of a two-particle 
state stems from the angular momentum: States with orbital angular momentum quantum number 
I involve the spherical harmonic Y^{n), where n = p/\p\ and p is the relative momentum of the 
two particles considered. Since y^(n) = (— l)^y^(— n), states with odd I have P and CP quantum 
numbers —1, while those with even / are even under P and CP. Since the decaying Kaon has no 
spin and the total angular momentum is conserved in any decay process, the two pions in the final 
state have have / = in the Kaon rest frame. (In general the spin wave function also matters, but 
pions have spin zero.) In total we find that the two-pion state is CP-even. Now ^Ciong was only seen 
to decay into three pions, so that this meson was believed to be CP-odd. In fact, its long lifetime 
stems from the kinematical suppression of the decay into the CP-odd three-pion state. To understand 
that a three-pion state is always CP-odd, first note that we get a contribution of (—1)^ = — 1 from 
the intrinsic CP quantum numbers of the three pions. Next pick any two of the pions and call there 
relative orbital angular momentum quantum number li. Likewise we denote the quantum number 
for the relative orbital angular momentum between this pair and the third pion by I2. One of the 
selection rules for the addition of angular momenta implies that the total quantum number / satisfies 
I > \li — hi- Since I = 0, this means that li = I2 and the "orbital" contribution to the CP quantum 
number is (— = (^_i)2«i _ ^ Thus the three-pion state is CP-odd, irrespective of the value of 
h. 

In 1964 the decay Kiong ^rvr was observed, establishing CP violation [4]. The two-generation 
Standard Model, whose construction was completed later in that decade [5], could not accommodate 
this phenomenon: We will see below that CP-violating interactions of quarks necessarily involve 
complex couplings. While the Vj^s in Eq. ^ are a priori complex, one can render them real in the 
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two-generation SM by transforming the quark fields as 

dj — > e'^idj, Uk —>■ e^'^i^Uk- (7) 

with appropriate phases (p'j and 0^. The net effects of these rephasings are the replacements of the 
Vjfc'sby 

^,.e^(^.-^^). (8) 

These expressions involves three independent phases and we may choose e.g. (pf — (p'f, (fif — (f)2 and 
^2 ~ 'Pi i'^ ^"'^h a way that the three complex phases of a unitarity 2x2 matrix are eliminated, arriving 
at the real Cabibbo matrix. In 1973 Kobayashi and Maskawa have pointed out that a physical CP- 
violating phase persists in the quark mixing matrix, if there are at least three generations [6] : A unitary 
3x3 matrix has 6 complex phases while we have only 5 phase differences (f)'^ — (f)^ at our disposal. 
The finding of Kobayashi and Maskawa was largely ignored at that time and only appreciated after the 
third fermion generation was experimentally established. In 1987 the ARGUS experiment at DESY 
observed B^ — Bd mixing, at an unexpectedly large rate [7]. This finding was the first hint at a truly 
heavy top quark, which enters the lower left box diagram of Fig. [T] 

1.3 CP violation 

The last stroke of the brush is devoted to CP violation. Defining 

CP\M[p^)) = -\M{-p^)), CP\M{p^)) = -|M(-%)) (9) 

we first look at decays M fcp andM fcp, where fcp is a CP eigenstate: 

CP\fcp) = ??cp|/cp> (10) 

with rjcp = ±1- The CP operator appearing in Eqs. ^ and (ITOl) is unitary, i.e. (CP)^^ = (CP)^ 
To get an idea of the importance of meson-antimeson mixing for the study of CP violation we first 
assume that AI and M do not mix. We could still measure the decay rates of the CP-conjugate 
processes M —>■ fcp andM — > fcp. If we find them different we establish direct CP violation 
(often called CP violation in decay). However, it is very difficult to relate a direct CP asymmetry to 
a fundamental CP phase in £: A non-zero direct CP asymmetry also requires final state interaction 
related to the rescattering process M f'^^ fcp- Rescattering leads to CP-conserving complex 
phases in the decay amplitude. In the absence of such phases the amplitudes of M — > fcp and 
M —>■ fcp would simply be related by complex conjugation since all phases would switch sign under 
CP. But then the two decay amplitudes would have the same magnitude leading to identical decay 
rates. For M = D^Bd, Bg this hadronic rescattering process is mainly inelastic and intractable with 
present theoretical methods. 

But thanks to meson-antimeson mixing we can study meson states which are superpositions of |M) 
and |M). The mass eigenstates \Mh) and \Ml) are linear combinations of |M) and [A/): 



\Ml) = p\M)+q\M), 
\Mh) = p\M)-q\M), 



(11) 



1.3 CP violation 
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with \p\^ + |gp = 1. We can calculate p and q from the box diagrams in Fig. [T]and will do so in the 
following sections. A commonly used shorthand notation for decay amplitudes is 

Af = A{M^f) = {f\S\M), Af = A{M^f) = {f\S\M). (12) 

A key quantity to study CP violation is the combination 

= (13) 

P 

Xf encodes the essential feature of the interference of the M ^ f and M ^ f decays, the relative 
phase between q/p (from meson-antimeson mixing) and Af/Af (stemming from the studied decay). 
In a first application, I discuss the decays of neutral Kaons into two charged or neutral pions. A neutral 
K ox K meson state is a superposition of Kh = -f^iong and Kl = ii'short- At short times the decays 
of the i^short component of our Kaon beam will vastly dominate over the Kiong decays and one can 
access the decay rates r(Ershort — > ^^r) for vrvr = vr+vr^, vr'^vr'^. At large times, say, after 10 times 
the -fCshort lifetime, our beam is practically a pure Kiong beam and we can study the CP-violating 
T{K\orig — > vrvr) decays. It is advantageous to switch to the eigenbasis of strong isospin /: 



TT TT ) = 



2 , . , , /T 



K^^ ) = Y 3 I (7rvr)^=o) + Y 3 I ^^^)^=2) ' 

The strong interaction respects strong-isospin symmetry to an accuracy of typically 2%, so that we can 
neglect any rescattering between the / = and 1 = 2 states. Consequently, no direct CP violation 
contributes to the famous CP-violating quantity 

_ ((^^)/=0|-^long) ^^4^ 
((vr7r)/=o|-f^short) 

Abbreviating Aq = Aq = and (see Eq. ([T3]l) Aq = \TnT)i^o insert Eq. 

into Eq. ([141) and readily find 

The experimental value [8] 

e^P = e*"^^ (2.23 ±0.01) X 10"^ with = (0.967 ± 0.001) ^ . (16) 

therefore allows us to determine Aq, which in our example is apparently close to 1. In our case with 
l^ol = l^ol (absence of direct CP violation) we have |Ao| = \q/p\- With Eq. ([T5] ) we find 



1 r 1 1 / 

ex ^ - [1 - Ao] 2 ^ ^ ~ 



-ilmAo (17) 



up to corrections of order ej^. Remarkably, from the real and imaginary part of ex we infer two 
CP- violating quantities: 



i) the deviation of from 1 and 
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ii) the deviation of Im Aq from 0. 

The first quantity is independent of the studied final state / and codifies CP violation in mixing. The 
second quantity, ImAj, measures CP violation in the interference of mixing and decay or, in short, 
mixing-induced CP violation in the decay M f. 

In the case of D,B^ or Bg mixing studies one tags the flavour at some time t = 0. The corresponding 
meson states are called \M{t)) and \M{t)) and satisfy \M{t = 0)) = |M) and \M{t = 0)) = |M). 
For t > these time-dependent states are calculable superpositions of \M) and \M) and by observing 
the time-dependence of M{t) — > / we can infer A/. The presently most prominent application of this 
method is the precise determination of Im Aj in the decay J/ipK^hort by the B factories BaBar 

and BELLE. Needless to say that we will discuss this important topic in detail below. 

While C, P, and T are violated in nature, the combination CPT is a good symmetry. This CPT 
theorem holds in any local Poincare-invariant quantum field theory [9]. It implies that particles and 
antiparticles have the same masses and total decay widths. When applied to our mixing problem 
characterised by S in Eq. ([6]) the CPT theorem enforces Sn = ^22- However, while the CPT theo- 
rem implies Tiot{M) = TtotiM), one still has different time-integrated total decay rates for tagged 
mesons, dtT tot {M{t)) 7^ dtT tot {M{t)). This quantity is sensitive to the "arrow of time" and 
the difference rtot{M{t)) — Ttot{M{t)) measures CP violation rather than CPT violation. Through- 
out my lectures I assume CPT invariance and therefore identify CP symmetry with T symmetry. Still, 
experiments have tested the CPT theorem by probing Sn = S22 in K — K mixing. We may specu- 
late that Poincare invariance and CPT symmetry are violated by the unknown dynamics of quantum 
gravity. If we are lucky the size of CPT violation scales linearly in the inverse Planck Mass Mpianck- 
Interestingly, today's accuracy of the CPT test Sn = S22 is roughly M/^ /Mpianck- 



2 Second lecture: Time evolution 
2.1 Time-dependent meson states 

In the Schrodinger picture, the time evolution of a quantum-mechanical state l-^) = l^',* = 0) is 
given by \ip, t) = U{t, 0)|^), with the unitary time-evolution operator U{t, 0). Consider first the case 
of a weakly-decaying charged meson (i.e. , or B^), which cannot mix with other states. The 
corresponding state at t = 0, |M+), will evolve into a superposition of all states allowed by energy- 
momentum conservation. This class of states consists of the original meson state \M^) and all final 
states I/) into which Af+ can decay. Defining 

\M+{t)) = |Af+)(A/+|Z^(i,0)|M+) (18) 

we can write 

U{tMM+) = \M+{t)) + Y,\f){fmMM+). 

f 

In order to find |M+(t)) we take a shortcut, by employing the exponential decay law to deduce 

|M+(t)) = e-**^«*e-r*/2|^+) (19) 



2.1 Time-dependent meson states 
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Figure 2: Left: generic self energy S of a charged meson. Right: M^ — M^ mixing amplitude S12. 



in the meson rest frame. The first term is the familiar time evolution factor of a stable state with energy 
E = Mm- The second factor involving the total width F is understood by considering the probability 
to find an undecayed meson at time t: 

\{M+\M+{t))f = e~^* 

Whenever I work in the Schrodinger picture I normalise the states as (Af+|Af+) = 1. Since Mm — 
iF/2 is independent of t, we can compute it using the familiar covariant formulation of quantum field 
theory. The optical theorem tells us that Mm and — F/2 are given by the real and imaginary parts of 
the self-energy S (depicted in the left diagram of Fig.O, where 

-ii2.W-P>^ = '^^^"^-'Ifl^'^'^T (20, 

zivim 

(To be precise, the diagram in Fig. |2] corresponds to 2Mm'^, so that S = Mm — iT /2 has mass 
dimension 1.) From Eq. ( [19] ) we find 

i^|M+(t)) = [MM-i^)\M+{t)). (21) 
This equation can be generalised to a two-state system describing neutral meson mixing: 

.d_(\M(t))\ ^ .22) 
'dt\\M{t))] \\M{t))] 

where now S is the 2x2 matrix defined in Eq. Recall that any matrix can be written as the sum 
of a hermitian and an antihermitian matrix. We write 

S = M-i- (23) 
2 

with the mass matrix M = and the decay matrix F = Then 

M,, = ^11±^, = (24) 



The expressions on the RHS of Eq. (1241) are called dispersive and absorptive parts of S12, respectively. 
The right diagram in Fig. |2] generically represents all contributions to S12. To compute E12 we 
can certainly use perturbation theory for the weak interaction (which to lowest order amounts to the 
calculation of the box diagram in Fig. [T]), but we must take into account the non-perturbative nature 
of the strong binding forces. The diagonal elements Mu and M22 are the masses of M and M and 
are generated from the quark mass terms in C and from the binding energy of the strong interaction. 
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However, the off-diagonal elements M12 = and all elements of T stem from the weak interaction 
and are therefore tiny in comparison with Mn and M22- The only reason why we can experimentally 
access M12 roots in the CFT theorem: CFT symmetry enforces 

Mil = M22, Til = r22, (25) 

so that the eigenvalues of S are exactly degenerate for S12 = S21 = 0. Even the smallest S12 can lift 
the degeneracy and can lead to large meson-antimeson mixing. 

With our shortcut we have avoided to prove that Eq. (|2TI ) holds with time-independent M and F. In 
fact, Eq. (|2TI) and the inferred exponential decay law in Eq. ([T9l ) are not valid exactly, but receive tiny 
(and phenomenologically irrelevant) corrections [10]. The same statement is true for Eqs. (l22l ) and 
(I23] ). a proper derivation of Eq. ((22l) using time-dependent perturbation theory for the weak interaction 
employs the so-called Wigner-Weisskopf approximation [11]. Corrections to this approximation have 
been addressed in Ref. [13] and are below the 10^^'^ level. 

We now proceed with the solution of our Schrodinger equation in Eq. (|22l) . Eq. (fTTI) means that the 
eigenvectors of S in Eq. (O are (p, q)^ and (p, —q)^. That is, S is diagonalised as 



with 



Q = ( ^ M and = — ( ^ M . (27) 

\Q -Q J 2pq\q -p J 

The ansatz in Eq. (|27] ) works because Sn = ^22- The mass eigenstates |Mi //(t)) obey an exponen- 
tial decay law as |M+(t)) in Eq. ([T9l ) with (Ma/ , F) replaced by {Ml^h,^l,h)- Transforming back 
to the flavour basis gives 

( ""(^^^ ] - q( ^""^"'^^'^ ' ]q~^( ^ (28) 

V \M{t)) J \ e-^M„t-^Ht/2 y ^ V 1^) / 

I adopt the following definitions for the average mass and width and the mass and width differences 
of the mass eigenstates: 

m = = Mil = M22 , r = = Til = r22 , ^29) 

AM = Mh-Ml, at = Tl-Th- 

Note that AM is positive by definition while AT can have either sign. Experimentally the sign of AT 
is only known for Kaons and my sign convention in Eq. ( [29l) corresponds to AFx > 0. The Standard- 
Model prediction for AF and AF is also positive, while no reliable prediction is possible for the 
sign of AF/). The matrix appearing in Eq. (l28l) can be compactly written as 




2.2 AM, Ar and CP violation in mixing 
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with 



9+it) 
9-{t) 



, Art AMt , Art AMt 
cosh cos I smh sm 



Art 

sinh — : — cos 



2 

AMt 



, Art 

+ i cosh — : — sm 



2 

AMt 



4 2 4 2 

Inserting Eq. (1301 ) into Eq. (1281 ) gives us a transparent picture of the meson-antimeson oscillations: 

|M(t)) 

|M(t)) 



(31) 



g+{t)\M) + ^g^{t)\M) 
P 



P 



^g_{t)\M)+ g+{t)\M) 



(32) 



We verify g+{0) = 1 and g-{0) = and find that g±{t) has no zeros for t > if Ar / 0. Hence 
an initially produced M will never turn into a pure M or back into a pure M. We will frequently 
encounter the combinations 



\9±{t)\' 

glit)g.{t) 



-rt 



Art 

cosh — ± cos (AM t) 



-rt 



sinh 



2 

Art 



+ z sin (AMt) 



(33) 



2.2 AM, Ar and CP violation in mixing 

We still need to solve our eigenvalue problem. The secular equation for the two eigenvalues ai^n = 
Ml,h — i^L,H/'^ of S is (Sii — (Jl,h)'^ — 5]i25]2i = 0. The two solutions of this equation therefore 
satisfy 

{aH-crhf' = 4i;i2S2i 

or 

Ar / r \ / r* \ 

(AM + i—f = 4 (^Mi2 - i^j [Mt, - i^j . (34) 
Taking real and imaginary part of this equation leads us to 

{AMf -^{AVf = 4|Mi2p - |ri2p , (35) 

AMAr = -4Re(Mi2n2)> (36) 
Further Eq. (l26l ) implies [Q^^TiQ]i2 = [Q^^^Q]2i = 0, which determines 

q ^ AM + iAr/2 _ 2Mf2-in2 

p 2Mi2-iTi2 AM + iAT/2' ^ ' 

(There is also a second solution with the opposite sign, which, however, is eliminated by imposing 
AM > 0.) For the simplification of Eqs. (l35l - [37l) it is useful to identify the physical quantities of the 
mixing problem in Eqs. (l22l) and (l23l) . In quantum mechanics we can always multiply either |M) or 
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|M) by an arbitrary phase factor without changing the physics. This will change the phases of M12, 
and q/p, none of which can therefore have any physical meaning. The three physical quantities 
of meson-antimeson mixing are 



12 



121 



and 



arg 



M12 



12 



Eq. (1361 ) then reads 



AM AT = 4|Mi2||ri2|cos0. 



(38) 



(39) 



We can easily solve Eqs. (1351 ) and ( [39l ) to express AAf and AT, which we want to measure by studying 
meson time evolutions, in terms of the theoretical quantities |Mi2|, |ri2| and (j). We recognise that the 
phase (f) is responsible for CP violation in mixing introduced after Eq. ([17] ): By multiplying the two 
expression for q/p in Eq. (|37] ) with each other we find 









l + i 


T12 










2Mi2 


2Mi2 




M12 


l + i 


ri2 


e-i0 








2Mi2 



(40) 



We immediately verify from this expression that cj) ^ Q,t: indeed implies \q/p\ / 1, which defines 
CP violation in mixing. 

Interestingly, CP violation in mixing is small (if quantified in terms of \q/p\ — 1) for the K, and 
Bg systems. For D — D mixing this is most likely also the case, but the experimental data are not 
accurate enough at present. In the case of K — K mixing we have established this phenomenon in 
Eq. (fTTl) from the measured value of Re ex in Eq. ([T6l ). In the B—B systems the line of arguments is 
as follows: Experimentally we know AM ^ AT and theoretically | ri2 1 <C AM is firmly established 
from a SM calculation, since the possible impact of new physics on |ri2| is small. Then Eqs. (1351 ) 
and ( [39l ) imply AM 2|Afi2| and therefore |ri2| ^ IM12I, so that the second term in the numerator 
and denominator of Eq. (l40l) is small, irrespective of the value of (p. Thus \q/p\ ~ 1 for Bd and Bs 
mesons. It is useful to define the quantity a through 



1 



(41) 



For the K, B^ and Bs systems we know that a is small. By expanding {q/p)'^ in Eq. (l40l ) in terms of 
(j) or ri2/Mi2 we find 



4|ri2| IM12I 

4|Mi2|2 + |ri2|2 

Infill + 



+ ( (Im " / 

M12 M12J 
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M12 



sm ( 



for _ftr-K mixing (42) 
for 5- ;B mixing. (43) 



With this result it is straightforward to solve Eqs. (1351 ) and (|39l ) for AM and AT. Incidentally, in both 
cases we have 



AM 
AT 



2IM12I, 
2 |ri2| cos ( 



(44) 
(45) 
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which holds up to corrections of order 0^ for Kaons and of order |ri2/Mi2p for B mesons. Of 
course, in the former case one can also replace cos by 1 . Importantly, in B physics one deduces 
from Eq. ^ that 



q _ Mi*2 



P \Mi2\ 



[1 + 0{a)] . (46) 



That is, the phase of —q/p is essentially given by the phase of the B^ — Bd or Bg — Bg box diagram 
in Fig.[T] Since B — B mixing is dominated by the box diagram with internal tops we readily infer 

^ = --i^ = -exp[iarg(y;,Ft,)2] for - mixing with g = d, s (47) 

P nbVtq 

up to tiny corrections of order a. 



2.3 Time-dependent decay rates 

Flavour factories are e+e" colliders whose CMS energy matches the mass of an excited quarkonium 
state which predominantly decays into (M, M) pairs. Running on the '0(3770), T(45') or T{5S) 
resonances, one copiously produces {D,D), {Bd,Bd) or {Bs,Bs) mesons. The {M,M) pairs are 
in an entangled quantum-mechanical state until the decay of one of the mesons is observed. If the 
decay mode M — > / is allowed while M — > / is forbidden one calls M ^ f a. flavour-specific 
mode or a tagging mode. The most prominent examples are the semileptonic decays M Xi^iy£. 
For the discovery of Bs—Bg mixing the flavour-specific mode Bs D^tt^ has played an important 
role [14]. A flavour-specific decay tags the decaying meson as either M or M. The Einstein-Podolsky- 
Rosen effect then ensures that the other meson is an M or M, respectively. The time of the flavour 
tagging "starts the clock", i.e. defines t = in Eqs. (|3T]) and (l32l ). This method is called opposite-side 
tagging. In hadron colliders pairs of different hadrons can be produced, e.g. a Bg can be produced 
together with a -B^ or A;, plus several lighter hadrons. Still, at the quark level (b, b) pairs are produced, 
so that the flavour tagging works as well. As an additional possibility, hadron colliders permit same- 
side tagging, where the flavour is determined at the time of the hadronisation process: When, say, a 
6-quark hadronises into a B meson several pions and Kaons are produced as well. The charges of 
these light mesons are correlated with the charge of the light valence quark, which in the case of the 
B meson is an anti-d quark. 

The time-dependent decay rate of a meson tagged at t = as M is defined as 

where dN{M{t) /) denotes the number of decays into the final state / occurring within the time 
interval between t and t + dt. Nm is the total number of M's produced at time t = 0. An analogous 
definition holds for T{M{t) /). One has 

r(M(f) ^ f)=Mf\{f\S\Mm\ r(M(t) ^f)=Mf \{f\S\Mit))\^ (49) 

with the time-independent normalisation factor Mj comprising the result of the phase-space integra- 
tion. It is straightforward to calculate r(M(t) f) and r{M{t) f) in terms of Aj and Aj 
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defined in Eq. we just need to insert \M{t)) and \M{t)) from Eq. ^ into Eq. (|49l). Trading Af 
for Xf (see Eq. ([T3] )) and a (see Eq. (|4TI )) and making use of Eq. (|33] ) we find the desired formulae: 



Art 1 - lA/i 



r(M(i) -^f) = Mf \A '2 




-Re A/ sinh— ^ him A/ sin(AMt) 



Often we want to compare these decay modes with the corresponding decays into the final state which 
is CP-conjugate with respect to /. For states / with two or more particles we define 

\7) = CP\f), (52) 

while for the initial one-particle states we have defined CP in Eq. Q. For example, for / = D^tt^ 
the CP-conjugate state is / = Df-K'. Whenever we discuss CP (or any other discrete transforma- 
tion) in decay processes, we apply the transformation in the rest frame of the decaying meson. The 
transformation in Eq. (|52l ) is understood to reverse the signs of three-momenta as in Eq. For two- 
body final states, which are our prime focus, we can rotate this mirror-reflected state by 180°, so that 
the three-momenta of the rotated CP-transformed state coincide with those of the original state. This 
procedure is usually implicitly understood when people discuss decays into CP eigenstates composed 
of two distinct particles, such as K ^ vr+vr". For a CP eigenstate fcp Eqs. (ITOl) and (l52l ) imply 
I7cp) = Vfcplfcp)- 



In the M{t) — > / decay rates it is advantageous to keep Aj while trading Aj for Aj 



— 2 fl + |A-F|-2 /^p^ l_|A-^|-2 

r(M(t) ^ /) = Mf Aj e"^* (1 - a) <( cosh — cos(AMt) 



-Re sinh + Im sin(AM t) > , (53) 

Aj 2 Aj J 

— - - 2 „ fl + IAjl^^ l-|Aj|-2 

r(M(t) ^ /) = Mf Aj e-^* <^ ^-i— cosh -— + cos(AMt) 

\ 2 2 2 

-Re sinh - Im sin( AM t) I . (54) 

Aj 2 Aj I 



Eqs. (|50] - [5TI ) and Eqs. (l53]454l ) are our master formulae to calculate any time-dependent decay rate 
of interest. We discuss two important applications here. The first one is the time dependence of a 
flavour-specific decay, which satisfies Af = Aj = Aj = 1/Aj = 0. In addition we consider a 
decay mode with | Aj\ = \Af\, that is without direct CP violation. Semileptonic decays satisfy both 
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conditions. Our master formulae become very simple for this case. Defining the mixing asymmetry, 

r(M(t)^/)-r(M(t)^/) (55, 
r(Af(t)^/) + r(Af(()^/) 



one finds to order a: 



_ cos(AMt) a 
^" cosh(Art/2) + 2 



^ cos2(AMt) 



(56) 



cosh^ (Art/2). 

Note that Ao{t) is not a CP asymmetry. Instead T{M{t) ^ /) oc |(M|M(t))|2 is proportional to 
the probabihty that an "unmixed" M decays to / at time t, while T{M{t) ^ /) cx \{M\M{t))\^ 
is the corresponding probability for the process M ^ M ^ f. The asymmetry Ao{t) is often 
employed to measure AM. In the ARGUS discovery of Bd—Bd mixing [7] no time-dependence was 
observed. Instead so-called like-sign dilepton events were observed in semileptonic {B^, B^) decays, 
meaning that one of the two mesons must have mixed. By counting these events and comparing the 
number with the number of opposite-sign dilepton events one can infer the quantity x = AM /T. The 
corresponding formula can be found by integrating our master formulae over t. 

The CP asymmetry in flavour-specific decays (often called semileptonic CP asymmetry) reads 

r(A/(t) ^ /) + r(M(t) ^ /) l + {l-af ' 
Define the untagged decay rate 

r[f,t] = T{M{t)^ f) + T{M{t)^ f) (58) 

to find: 

_ r[f,t]-T[f,t] _ gfe gfe cos(AMt) 
afs,unt(tj - - 2 2 cosh(Art/2)- ^^^^ 

Hence no tagging is needed to measure Ofs! We observe that we can determine the three physical 
quantities characterising meson-antimeson mixing, |Afi2|, |ri2| and a, by measuring AM, AT and 
Ofs. At present all three quantities are only measured for K — K mixing! Also the semileptonic CP 
asymmetry of B mesons can be measured without observing any time dependence. In the spirit of 
ARGUS we can compare the number of positively-charged like-sign dilepton pairs with the number 
of negatively-charged ones. Such measurements are performed at the B factories and the Tevatron, 
but no non-zero semileptonic CP asymmetry has been established by now. 



Amusingly, the oscillations drop out from the tagged quantity in Eq. (1571) . while they persist in 
Eq. ( [59l ). In most applications one can neglect the tiny a in Eqs. (l50l - [5TI ) and Eqs. (l53]454l) . Then 
we realise that in the untagged rates, obtained by adding Eqs. (l50l) and (|5T]) or Eqs. (l53l) and (l54l ). the 
terms involving cos(AMt) and sin(AMt) vanish. 

The second application of our master formulae are decays into CP eigenstates, M — > /cp. The 
time-dependent CP asymmetry is 

... _ r(M(t)^/cp)-r(M(t)^/cp) 

- T{M{t) - /cp) + r(M(t) ^ /cp) • ^'"^ 



16 



3 Third lecture: Linking quarks to mesons 



Using Eq. (l50l) and Eq. (|5TI ) one finds 

^g^cos(AMt) + Aggsin(AMt) 
- -cosh(Art/2)+^Arsinh(Art/2) + ^^^^ ' ^^^^ 

with (for / = /cp) 

2 

^jj. 1 — I A J I jj^j^ 2ImAj 2ReAj 

" 1 ^ IX |2 ' ^CP = ,2 ' ^Ar = - I |2 - (62) 

1 + IA/l 1 + 1 + 1-^/1 

Note that \A'^''p\'^ + I^^J^P + \A/^t? = 1- Experimentally one can track the time-dependence of 
af{t) and read off the coefficients of cos(AMt) and sin(AMt), so that one can determine |A/| and 
Im Aj. When studying decay amplitudes we can treat the weak interaction perturbatively by drawing 
quark-level Feynman diagrams involving the exchange of W-bosons. While we cannot fully compute 
those diagrams, because we cannot estimate how the quarks are "dressed" by the strong interaction, we 
can still assess the CP-violating phases by identifying the CKM elements in the diagrams. Decays in 
which all contributing Feynman diagrams carry the same CP- violating phase are called golden modes. 
These modes satisfy \Af \ = \ Aj\ , so that there is no direct CP violation. In a golden M fcp decay 
this means | A \ = 1 and in Eqs. (|6T]) and (l62l) we have A^p = and 

^gj5 = ImA,fcp. (63) 

Moreover the phase of Afi^p/Afi-.p is trivially read off from the phase of the CKM elements. In B 
physics, where we also know the phase of q/p from Eq. (|47] ). we can therefore directly relate the 
measured ImXj^p to phases of CKM elements, if M ^ /cp is golden. 



3 Third lecture: Linking quarks to mesons 
3.1 The Cabibbo-Kobayashi-Maskawa matrix 

We have encountered the CKM matrix V in Eq. ([3]). A unitary 3x3 matrix can be parameterised by 
three angles and six complex phases. With the rephasings in Eqs. Q and ^ we can eliminate five 
phases from V leaving us with one physical CP-violating phase. In the parameterisation favoured by 
the Particle Data Book one has 

C12C13 S12C13 
V = I -S12C23 - ci2S23Si3e*'^i3 C12C23 - si2S23Si3e*'^'^ S23C13 | , (64) 

S12S23 - Cl2C23Sl3e*'^'^ -C12S23 - Sl2C23Sl3e*'^" 

where Cij = cos 9ij and Sij = sin Oij. The real angles 9ij may be chosen so that < 9ij < 7r/2, and 
the phase ^13 so that — vr < ^13 < vr. For the discussion of CKM metrology it is useful to introduce 
the Wolfenstein parameterisation [15] 

l-iA2 A AX^{p-ir]) 

V = \ -A 1 - AA2 I + O(A^) , (65) 

AX^{1- p-irj) -AX'^ 1 




3.1 The Cabibbo-Kobayashi-Maskawa matiix 
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which is an expansion in terms of the small parameter A = 0.22. The remaining three parameters 
A, p and are a bit smaller than 1. The Wolfenstein parameterisation nicely reveals the hierarchical 
structure of the CKM matrix, with diagonal elements of order 1 and smallest elements in the upper 
right and lower left comers. We can now understand why the prediction of rric from IS.Mk in 1974 was 
successful: Any contribution involving the top quark (at that time unknown and unimagined by the 
authors of Ref. [3]) to the upper left diagram in Fig. [His highly suppressed by small CKM elements, 
since |VtdVjs| ~ A^, while |VcdVcs| — — sin 0c — A. Further the upper left 2x2 submatrix, 

the Cabibbo matrix, is almost unitary and involves only a single parameter, the Cabibbo angle 9c with 
Vud — Vcs — cos^c and Vus — —Vcd — Therefore the two new elements Vd and Vcs predicted in 
Ref. [2] were completely fixed in terms of the known Oc- In the Wolfenstein approximation only Vub 
and Vtd have a complex phase and CP violation is characterised by r/ / 0. 



Any unitary 3x3 matrix satisfies 



and 



v;,Vki + v*^Vk2 + v*3Vk3 = Sjk- 



(66) 
(67) 



If we choose j ^ k the three terms add to zero. We can depict the relations in Eqs. (1661 ) and (1671 ) 
as triangles in the complex plane, e.g. for Eq. (l66l) the three comers are located at 0, V^jVik and 
— ^2* V2fe. The three sides can be associated with the three terms summing to zero. The area of all six 
triangles is the same and given by J/2, where J is the Jarlskog invariant [16] 



J 



Im [V:dVtbV:,V, 



ud\ 



C12C23C13 



Si2S23Sl3Sin(5i3 



(68) 



Here the third expression refers to the exact parameterisation of Eq. (l64l ) and the last result uses the 
Wolfenstein approximation. Four of the six unitarity triangles are squashed, the thi^ee sides are similar 
only for the choice {j,k) = (3, 1). Moreover, within the Wolfenstein approximation the shapes of 
the triangles corresponding to Eqs. (l66l ) and (l67l) are equal for (j, fc) = (3, 1). Applying the phase 
transformations of Eqs. (O and ([8]) rotates the unitarity triangles in the complex plane, but leaves their 
shape fixed. Seeking a definition of a rephasing-invariant unitarity triangle with a physical meaning 
we divide Eq. ([66]) (for (j, k) = (3, 1)) by 1^23 ^21 = V*^Vcd to arrive at 



v:,Vcd Y:,Vcd 







(69) 



When people speak of "the" unitarity triangle they refer to the rescaled triangle defined by Eq. (1691 ). 
Since its baseline coincides with the interval [0, 1] of the real axis, the unitarity triangle is completely 
determined by the location of its apex (p, r/), where 



p + iri 



VlVcd 



(70) 



Inserting Eq. (l65l) into Eq. (iTOl ) one realises that (jj, r]) = (p, rj) within the Wolfenstein approximation, 
which here is good to an accuracy of 3%. The unitarity triangle is depicted in Fig. [3] The two non- 
trivial sides of the triangle are 



CP-violating quantities are associated with the triangle's three angles 



a 



arg 



VudV. 



13 = arg 



uh - 



VtdV, 



7 = arg 



tb- 



^udVub 

VcdV^b. 



(71) 



(72) 



18 



3 Third lecture: Linking quarks to mesons 



p-nri 




l-p-ir| 



C=(0,0) B=(1,0) 
Figure 3: The (standard) unitarity triangle. 



The angle 7 coincides with ^13 of Eq. (I64b at the sub-permille level. With Eqs. 



p + irj = Rue^'^, I - p - irj 

The unitarity relation of Eq. ( [69l ) now simply reads 



Rte 



one obtains 
(73) 

(74) 

Taking real and imaginary parts of Eq. (1741) reproduces formulae which you know from high-school 
geometry, allowing us to express any two of the four quantities ii^, 7, /3 in terms of the remaining 
two ones. By multiplying Eq. d74l ) with either exp(— ^7) or exp(i/3) one finds analogous relations 
involving a = vr — /? — 7. 

Sometimes one needs to refine the Wolfenstein approximation to higher orders in A. It is prudent to 
define [17] 



A = S12, 



AX' 



S23 



(75) 



to all orders in A and to expand all CKM elements in terms of A, A, p and to the desired order in A. 
Then, for example: 



= AX^-p-irDil + ^ + OiX") 



The phase 



(is = arg 



A277 + 0(A^) 



(76) 



(77) 



plays an important role in Bg—Bs mixing; Ps is small, of order 0.02 (equal to 1 degree). In the phase 
convention of Eq. ( |64l ) the phase of VcsV^ is O(A^) and 



a.rgi-Vts) = Ps{l + 0{X^)). 



(78) 



Organising the phases in powers of A, we find all CKM elements real to order A^ except for Vub, Vtd 
and Vts- Going to higher orders one encounters arg(— V"cd) ^ A'^fjX^ and arg(ycs) — —A'^fjX^. 

3.2 Effective Hamiltonians 



We now address the strong interaction, which is the main obstacle on our way from quark diagrams 
to mesonic amplitudes like M12 and A{M — > /). In Sect. II. II we have seen that weak processes of 



3.2 Effective Hamiltonians 



19 




q 

Figure 4: The four-quark operator Q for Bq — Bq mixing with q = dot s. 



mesons are multi-scale processes. For instance, B—B mixing involves three largely separated scales, 
since mt ~ Mv/ » S> Aqcd- These scales must be disentangled to separate the short-distance 
QCD, which is described by the exchange of quarks and gluons, from the long-distance hadronic 
physics, whose characteristic property is the confinement of quarks into hadrons. The key tool to 
separate the physics associated with the scale mheavy from the dynamics associated with mught 
mheavy is the construction of an effective field theory. The corresponding effective Hamiltonian H^^ 
is designed to reproduce the S -matrix elements of the Standard Model up to corrections of order 
("T-iight/^^heavy)" whcrc n is a positivc integer: 



/ flight \ 

\ ^hcavy / 



(79) 



I exemplify the method with an effective Hamiltonian which reproduces the amplitude for B — B 
mixing up to corrections of order m^/M^. That is, we employ Eq. ( 1791 ) for the case i = B and 
f = B (where B = B^or Bg), mught = rrib and mheavy = Mw ~ nit. The corresponding effective 
Hamiltonian reads 



H' 



eff 



jjQCD(f=5) _^ jjQED{f=5) _^ jj\AB\-- 



(80) 



Here the first two terms are the usual QCD and QED interaction Hamiltonians with 5 "active flavours", 
meaning that they do not involve the top quark. The last term describes the weak interaction. Adapted 
to the process under study, //I'^^N^ only encodes the physics related to B — B mixing, but does 
not describe other weak processes such as meson decays. It is called i^l^^N^^ because it describes 
physical processes in which the bottom quantum number B changes by two units. //I^^l=2 does not 
contain W-boson, Z-boson or top-quark fields, instead the AB = 2 transition of the box diagram in 
Fig. [T]is mediated by an effective four-quark coupling: 



Q 



with q = dox s. 



(81) 



For historical reasons Q is called & four-quark operator, but it is nothing but a point-like coupling of 
four quark fields as shown in Fig. |4] We have 



4^ 



(82) 



where the lengthy expression multiplying Q is just the effective coupling constant multiplying the 
four-quark interaction of Fig.|4] This coupling constant is split into several factors, the first of which 
contains the Fermi constant Gp- The second factor summarises the CKM elements of the box diagram 
and the third factor Cl'^^l^^(mt, Mw,lj) is the Wilson coefficient, which contains the information on 
the heavy mass scales Mw and mt. Finally is the renormalisation scale, familiar from QCD. Just as 
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any other coupling also Q must be renormalised. The renormalised operator Q depends on /x through 
the renormalisation constant Zq{^) via Q = ZqQ^^^^ and (in a mass-independent scheme like MS) 
the latter dependence is only implicit through ^(/u), where g is the QCD coupling constant]^ With the 
decomposition in Eq. ( [82l ) C*!^^!^^ has dimension two and is real. 

(j\AB\=2 calculated from the defining property of H'^^ in Eq. (1791 ): We compute the AS = 2 
process both in the Standard Model and with the interactions of H^^ and adjust C''^^'^^ such that 
the two results are the same, up to corrections of order m^/M^. Obviously we cannot do this with 
mesons as external states i and /. But a crucial property of H^^ is the independence of the Wilson 
coefficient on the external states. We can compute it for an arbitrary momentum configuration for the 
external quarks as long as the external momenta are of the order of miight- That is, we do not need to 
know the complicated momentum configuration of quarks bound in a meson state. Further all QCD 
effects in C*!^^'^^ are purely perturbative: 



AB\ 



=2,(0) _^ 

4tt 



^|AiJ|=2,(l) _^ _ ^ _ 



(83) 



We can understand why and how this works if we expand the result of the box diagram of Fig. [T] in 
terms of the momenta of the external quarks, which are at most of order rrii,. The leading term consists 
of the result of a loop integral with external momenta set to zero and the spinors of the external quark 
states. Now the "effective theory side" of Eq. ( 1791 ) involves the tree-level diagram corresponding to 



AB1=2,{0) 



{fm: 



,{0) 



where \i) = \ph, Sfe;Pg, Sq) and |/) = \pg, Sq;p^, s^) are the external states characterised by the mo- 
menta and spins of the quarks. The superscript "(0)" indicates the lowest order of QCD everywhere. 
Since {f\Q\i) reproduces the spinor structure ("Dirac algebra") of the box diagram, the coefficient 
(j\AB\=2,{o) inferred from this matching calculation is solely determined in terms of the loop integral 
and therefore only depends on Mw and rrit. The matching calculation becomes less trivial when we 
go to the next-to-leading order (NLO) of QCD. Now i^Q^D enters the matching calculation and we 
must dress both the box diagram and the effective diagram in Fig. |4] with gluons in all possible ways. 
Denoting the SM amplitude by 



M 



4tt 



(84) 



our NLO matching calculation amounts to the determination of C'^^' ^'^^^ from 



M 



(0) 



An 



47r2 



tqj 



q\AB\=2,(0) _^ ^q\AB\=2,(1) 



(g)(0) + 



l + O 



+ O (a2) (85) 



On the RHS the external states are suppressed for simplicity of notation. The QCD corrections to 
the box diagram in TW^^^ not only depend on the light scales, i.e. external momenta and light quark 



*The analogy with the renormalisation of the QCD coupling constant is more obvious if one reads the product CZqQ'""^^ 
in a different way: By assigning Zq to C rather than Q one may view C as a renormalised coupling constant. The notion 
of a "renormalised" operator instead of a "renormalised Wilson coefficient" has historical reasons. 



3.2 Effective Hamiltonians 



21 



masses, they also suffer from infrared (IR) divergences. These divergences signal the breakdown 
of QCD perturbation theory at low energies. However, the gluonic corrections to Fig. |4l which are 
comprised in {Q)^^\ exactly reproduce the infrared structure of the SM diagrams: They involve the 
same IR divergences and have the same dependence on the light mass scales. Collecting the 0{as) 
terms from Eq. 



■2 



47r^ ^ L 



(86) 



one finds identical IR structures on the LHS and in the first term in the square brackets, while 
(j\AB\=2,{i) Qj^jy contains heavy masses and no IR divergences. In conclusion, the IR structure of the 
SM amplitude properly factorises with an "infrared-safe" C'l'^^l^^. This success can be understood 
by separately discussing the regions of small and large loop momentum passing through a gluon line 
in the diagrams of The infrared-sensitive diagrams are identified as those in which the gluon 

connects two external quark lines. (The other diagrams are infrared-finite and one can set the light 
mass parameters to zero.) If the loop momentum traversing the gluon line is small, we can neglect it in 
the heavy top and W propagators. Therefore the loop integration factorises into two one-loop integra- 
tions and the second loop integral involving the heavy particles simply reproduces the one-loop result 
contained in C^^^^^'^'^^\ The gluon-loop integration — still over soft momenta only — is equal to the 
one in the corresponding diagram in {Q)^^\ where the gluon connects the same quark lines. Therefore 
the region of integration with a soft gluon factorises with the leading-order coefficient C^^^\='^'('^) in 
Eq. (1851) . The region of the momentum integration with a hard gluon momentum does not factorise 
in this way and contributes to C^^^l^'^'W _ However, the region of large gluon loop momentum is 
not infrared-sensitive and we can neglect the light momenta and masses. Therefore c\^b\=2,{i) ^^^^ 
not depend on the light mass scales. Conversely, (Q) contains only small scales of order mught and 
encodes the full infrared structure of A4. Therefore our quark-level calculation is meaningful for 
(j\AB\=2^ but not for {Q). In order to make a theoretical prediction for the B — B mixing amplitude, 
we must compute {B\Q\B) with nonperturbative methods. The factorisation of M. into short-distance 
coefficients and long-distance operator matrix elements is also called operator product expansion. 

Here I only derive the result for the leading-order (LO) Wilson coefficient C'l'^^l^^'*^*^-'. In a first step 
let us decompose M '•^^ as 

j,k=u,c,t 

where -M^j^j^ {Q)^^^ is the result of the box diagram containing internal quark flavours {j, k) with the 
CKM elements factored out. We then write 

Mf^ = -^M^S{x„Xk) (88) 

with Xj = mj/My^^. The function S{xj,Xk) is symmetric, S{xj,Xf:) = S{xk,Xj). In the next step 
we use CKM unitarity to ehminate V*^Vuq = -Vtiytq - V*^Vcq from Eq. 

_A^(o) = ^Mir \{y:,vtqf s{xt) + 2v:,vtqv:,v,qS{x,,xt) + {v:,Vcqf s{xS\ 



47r2 

S and S are related as 



(89) 



S{xj,Xk) = S{xj,Xk) — S{xj,0) — S{0,Xk) + S{0,0), for j,k = c,t, 

S{x) = S{x,x), (90) 
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where I have set the up-quark mass to zero. In Eq. (|89l ) the last two terms are tiny, because Xc ~ 10 ^ 
and 

S{x,) = 0{xc), S{x„xt) = 0{x,lnx,). (91) 

This consequence of CKM unitarity is called the Glashow-Iliopoulos-Maiani ( GIM) suppression, re- 
lated to the vanishing of FCNCs in the limit of equal internal quark masses (here rric and m„ = 0). 
No GIM suppression occurs in top loops, because xt ~ 4. The dominant contribution to Eq. (l87l) 
involves 



S{xt) = Xt 



19 1 3 1 

+ 



4 Al-xt 2(l-xt 



\2 



1 3 



l-Xt 



Inxt ^ 2.3. (92) 



The tiny charm contribution does not contribute to C'l^^l"^'^'^) at all; to accommodate for it we must 
refine our operator product expansion to include higher powers of (w-ught/^T^heavy) in Eq. (1791 ). We 
can read off cl^-B|=2,(o) fj-om Eqs. ^ and §9^: 



Cl^^l=2'(0)(r„t,iV%,;u) = M^S{xt). (93) 
The functions S{x) and S{xc, xt) are called Inami-Lim functions [28]. 



The factorisation in Eqs. (1791 ) and (1851 ) also solves another problem: No largely separated scales 
appear in (^'^^'^^(mt, Mv/,n) provided that we take ^ = 0{Mw ,mt), so that no large logarithms 
can spoil the convergence of the perturbative series. While no explicit ^u-dependence is present in our 
LO result in Eq. (l93l ). there is an implicit ^u-dependence through mt(/i), which is a running quark 
mass (typically defined in the MS scheme). C^'^^\='^'W also contains an explicit ln{^/Mw) term. 
Two sources contribute to this term: First, there is already a ln{n/Mw) term in ^A^^\ familiar to 
us from matrix elements with MS-renormalised UV divergences. Second, M^^^ contains the large 
logarithm ln{mi,/Mw) which is split between matrix elements and Wilson coefficients as 

In^ = In^+ln-^. (94) 



This feature is transparent from Eq. (I86l) . 

The scale i^tw = O^My/ ^mt) at which we invoke Eq. (l85l) to find C'l'^^'^^ is called the matching 
scale and Cl^^l^^(m(, Mvi/, ^itw) has a good perturbative behaviour. Similarly, no large logarithms 
occur in {Q{^b)), if we choose a scale fn, ~ rrih in the matrix element. Since the ^u-dependence in 
jj\AB\=2 spurious, we can take any value of fi we want, but this value must be the same in C(/i) and 
{Q{^)). That forces us to either relate C{fitw) to or to express {Q{fii,)) in terms of {Qi^J-tw)) 

in such a way that large logarithms 

a^ln"^ (95) 

are summed to all orders n = 0,1,2... in perturbation theory. This can be achieved by solving 
the renormalisation group (RG) equation for either C(^) or {Q{ii)). All steps of this procedure are 
analogous to the calculation of the running quark mass, which can be found in any textbook on QCD. 
RG-improvement promotes our LO result to a leading-log (LL) quantity: 

(Qil^tw)) = n^^Hl^b,l^tw){Q{lib)) (97) 
u'-'''{fJ-b,fJ'tw) = I ^ ^^^^ J ° With 7^=4. (98) 



3.2 Effective Hamiltonians 
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The evolution factor u^^\fj,b, fxtw) depends on the anomalous dimension of Q, which equals 

^0 



(aj(47r))7|°^ to LL accuracy. P^/^ = 11 - 2//3 is the first term of the QCD /3 function. One 



usually writes 

C\^^\=\mt,Mw,fJ^b) = VBbB{f^b)C\^^^=^'^'^\rnt,Mw,l^tw) (99) 

where all dependence on /x;, is absorbed into bsi^J-b) and all heavy scales reside in rjs- This factorisa- 
tion is possible to all orders in a^. It is trivially verified in the LL approximation of Eq. ( |98l ). where 
simply m(°)(//6, //tvy) = VBbBifJ'b)- In Eq. ( |99l ) mt is understood as ■mt{mt) (and not as rtitintw))- 
In this way r^s is independent of ^tw to the calculated order; the residual fitw dependence is already 
tiny in the NLL result, tjb mildly depends on xt = rn^/My^ and in practice one can treat it as a 
constant number [18]: 

rjB = 0.55, hsiHb = mb = 4.2 GeV) = 1.5. (100) 

The dependences of bs on fib and the chosen renormalisation scheme cancel in the product 
bBif'b){Q{fJ'b))- The quoted number is for the MS-NDR scheme, where "NDR" refers to the treatment 
of the Dirac matrix 75. Details on this topic can be found in [19]. We see that the impact of short- 
distance QCD corrections is moderate, since r]B bBifJ-b) = 0.84. The NLL calculation of Ref. [18] has 
found only small two-loop corrections and the remaining uncertainty affects rjB only in the third digit 
behind the decimal point. RG-improved perturbation theory works superbly! Combining Eqs. (l82l ). 
( |93] ) and ( |99l ) we obtain our final expression for the | AB \ = 2 hamiltonian: 

^\AB\=2 ^ ^M^{VtbV*frjBS{xt)bB{fib)Q{f^b) + h.c. (101) 

Finally we cannot escape from quark confinement! Our hadronic matrix element is conventionally 
parameterised as 

{Bg\Q{,,b)\B,) = ^m2^4^^ (102) 

with the Bq meson decay constant fs, and the bag factor Bb,- The parameterisation in Eq. (11021) 
is chosen in such a way that BBq/bsilJ'b) is close to one. It will be especially useful once precise 
experimental data on /^^ ~ fB+ from leptonic decays will be available. With the help of our 
effective field theory we have beaten the problem of long-distance QCD mB—B mixing down to the 
calculation of a single number. Lattice gauge theory computations cover the ranges [29] 



fB,\/BB, = (225 ±35)MeV, fB^^/Bs, = (270±45)MeV. (103) 

The quoted hadronic uncertainties are the main problem in the extraction of |14{,14q| from the measured 
AMb,- Bb^ could differ from Bb^, but no computation has established any significant difference by 
now. 

Putting Eqs. (IIOII ) and (11021 ) together we find the desired element of the B — B mass matrix: 

{Bq\H\^^\=^\Bg) 



M12 



2Mb, 

r?B Mb, BbJI, ^^h/ ^(^) [VtbVt^f . (104) 



127r 
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Figure 5: Second-order contribution of ifl^^l ^ to Bg — Bs mixing. The diagrams constitute the 
dominant contribution to AT ■ 



We remark that there is no contribution of //l'^^l=2 ri2, because {Bq\H^'^^^^'^\Bq) has no absorp- 
tive part. By inspecting Eq. (l24l) we can verify that the dispersive or absorptive part of some amplitude 
can be calculated by replacing the loop integrals by their real or imaginary parts, respectively, while 
keeping all complex CKM elements. But only diagrams with light internal quarks involve loop inte- 
grals with a non-zero imaginary part. Hence we must extend our effective-Hamiltonian formalism to 
include the effects of light internal quarks in the box diagrams, if we want to predict Ar^^. Contract- 
ing the heavy W-boson lines in the diagrams of Fig.[T]to a point does not correspond to a contribution 
from i7l'^^l=2 jj^ tj^e effective theory. Instead this is a second-order effect involving some effective 
\AB\ = 1-Hamiltonian i/l^^l=i, which we must add to H^^ in Eq. (l80l ). The relevant piece from the 
RHS ofEq. (I79l)is 

- ^ y d^xd^y {B\TH^^^\=\x)H\^^\=\y)\B). (105) 

The LO contribution to this bilocal matrix element is depicted in Fig. [5] for the case of Bg — Bg 
mixing. The contribution from Eq. (11051 ) to B — B mixing is much smaller than the one from 
^|AB|=2^ which is enhanced due to the heavy top mass entering Eq. (l92l ). Therefore we can neglect 
the bilocal contribution in M12 and only need to consider it for ri2. From this observation we also 
conclude that |ri2| < IM12I leading to |Ar| < AM, which we akeady exploited in Eqs. (143^47] ). 



3.3 SM predictions of AM, AT and 

In Sec. 13. 21 we have collected all ingredients of the SM calculation of AM = 2|Mi2| for the B^ and 
Bs systems. Looking at Eq. (1651 ) we realise that | Vtb \ is well-known and | Vts \ is essentially fixed by 
the well-measured |V"cfe|. From Eqs. (11041) and (11031) we find the SM prediction 

AMb, = (12.5 ± 4.3) meV = (19.0 ± 6.6) ps"^ (106) 

The first unit is milli-electronvolt, a unit which we do not encounter often in high-energy physics. By 
dividing with h one finds the second expression in terms of inverse picoseconds, which is more useful 
since AM is measured from the oscillation frequency in Eq. (1551 ). Eq. (11061 ) is in good agreement 
with the Tevatron measurement of [14, 20] 

AM^^^P = (17.77 ± O.lO(stat) ± 0.07(syst)) ps"^ (107) 

The corresponding quantity for Bd — Bd mixing is well-measured by several experiments with [8] 

AM^^'P = (333.7 ± 3.3) //eV = (0.507 ± 0.005) ps"^ (108) 
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We can use AMb^ to determine | Vtd \ ■ From Eq. (11041) we infer 



AMb^ = (0.52 ±0.02) ps' 



\Vtd\ 
0.0082 



fBi\jBB^ 

225 MeV 



(109) 



The 16% error of the lattice value in Eq. (11031 ) dominates the uncertainty on the extracted \ Vtd\- The 
all-order Wolfenstein parameterisation defined by Eqs. (TTOI ) and (1751 ) implies 



(110) 



\Vtd\=A\'Rt + 0{\''). 



Since ~ |Fcfc| is well-known, AMb^ essentially determines Rt, i.e. one side of the unitarity 
triangle. Even better, we can use the ratio AM^^/AMs^ for the same purpose: If one forms the ratio 
of the hadronic quantities in Eq. (11031 ). many uncertainties drop out: 



fSs^BB 



1.20 ±0.06. 



(Ill) 



fSayBBa 



In the limit of exact flavour-SU(3) symmetry (corresponding to = = m^) one has ^ = 1 
which reduces the calculational task to compute the deviation of ^ from 1. The somewhat large error 
in Eq. (|llll) reflects the ongoing discussion on potentially large chiral logarithms [21] which may 
increase ^ significantly. This problem occurs, because lattice simulations use values for the pion mass 
which are larger than the physical value. The extrapolation to tUt^ ~ 140 MeV with the help of chiral 
perturbation theory introduces this source of error. Sum-rule calculations of (or rather fss/ fB^) 
which automatically include these logarithms, however, give values at the lower end of the range in 
Eq. (|llll) [22]. Further all short-distance QCD drops out from the ratio AMb^/ AMb^, so that one 
simply has 



Vtd 



Vt 



ts 



IAMb, Mb^ 



AMb, V Mb, 



(112) 



The Wolfenstein expansion leads to 

Vtd 



Vt 



ts 



RtX 



(113) 



Combining Eqs. (|112l) and (1113b (and using Mb^ I^b^ = 1-017) we easily derive a home-use formula 
for Rt. 



Rt 



0.887 



AMb, 17.77ps-^ e 



A 



0.507 ps-i AMb, 1.2 0.2246 



1 ±0.05p] 



(114) 



Neither p k, 0.2 nor the 1% error on A ~ 0.2246 have an impact on the error of Rt- Using the 
numerical input from Eqs. (1107141081 ) and Eq. (Illll ) we find 



Rt 



0.90 ± 0.04 



(115) 



and the uncertainty is essentially solely from ^ in Eq. (Illll ). 
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Next we discuss AF and the quantity Ofs in Eq. (|57] ). which governs CP violation in mixing. In order 
to find these quantities we need to calculate ri2. This involves the diagrams of Fig. [5] and brings in a 
new feature, power corrections of order AqcD/^Tife [23]. NLL QCD corrections to F12 in the B system 
have been calculated in Ref. [24-26]. In the SM the CP phase (j) of Eq. ( [38] ) is so small that one can 
set cos (/) to 1 in Eq. ([45] ). If we normalise AF to AM we can eliminate the bulk of the hadronic 
uncertainties. Updated values, obtained by using an improved operator basis, are [27] 

0.088 ±0.017ps-\ (116) 

(2Q.7tlfj-W~^ps-\ (117) 

The width difference in the Bs system amounts to 12.7 ± 2.4% of the average width Tb^ —^b^ [27] 
and is in the reach of present experiments. Needless to say that there are no useful data on AF^^. The 
predictions for the CP asymmetries in flavour-specific decays of Eq. (l57l) are calculated from Eq. (l43l) 
and read [25-27] 

a| = (2.06 ± 0.57) • 10"^ (118) 
ai = (-4.8^1:°) -10-^ (119) 

Also the current data for these CP asymmetries are not useful for CKM metrology. A future measure- 
ment of Ofg'^^^ will add an interesting new constraint to the (p, fy ) plane [25] : 

dexp 

(?7 - Ri,f + (1 - pf = Rl with = — ^p- . (120) 

(lO.lLnj ■ 10-4 

The theory prediction of Refs. [25,26] enters the denominator of Rf^, the quoted value is consistent 
with Eq. (II 191 ) and stems from the update in Ref. [27]. Eq. ( 11201 ) defines a circle with radius Rf^ 
centred around {p,f]) = (1, i?fs). Therefore the circle touches the p axis at the point (1,0), see Fig.|6] 

We have seen that the three quantities related to Bg — Bg mixing discussed in Eqs. (11061 ). (11161 ) and 
(11181 ) have little dependence on p and fj. Only AMb^ has an impact on CKM metrology, through 
Eq. (11141 ). The small sensitivity to p and fj becomes a virtue in searches for new physics, where 
Bg — Bg mixing plays an important role. 

Next we discuss K — K mixing: The calculation of M12 now forces us to compute box diagrams 
of Fig. [T] with all possible quark flavours u, c, t, because the top contribution involving S{xt) is sup- 
pressed by the small CKM factor {VigVtd)^ — A'^\^^{\ — p + vq)'^- The charm and up contributions, 
however, are proportional to only two powers of A. Therefore we cannot neglect these contributions 
despite of the smallness of S{xc) and S{xc-,xt) (discussed around Eq. (|9T])). Their calculation pro- 
ceeds in two major steps: First, the top quark and W-boson are integrated out. In the resulting effec- 
tive theory the AS" = 2 transitions receive second-order contributions from a |A5| = 1-Hamiltonian 
^|AS'|=i have, already seen this in our discussion of Ai? = 2 transitions, the corresponding 
expression iov K — K mixing is obtained by replacing ffl'^^l=^ with ffl^'^l=i in Eq. (11051 ) (and is 
described by the analogous diagrams of Fig.[5l). In addition to this bilocal contribution, the term with 
5(xc, xt) also involves a | A5| = 2-Hamiltonian H^^^^^'^ which mediates K — K mixing via a local 



AFb = ( AMS^P = 
AFb, = ( AM^^P = 
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Figure 6: Impact of af^ on the {p, fj) plane: The solid blue curves limit the allowed range (defined by 
the error in Eq. (11201) ) for a hypothetical measurement of a^^^^^ = —5 • 10~^. The solid red curves are 
for Ojg"^^^ = —10"'^ instead. For further information see Ref. [25], from which the figure is taken. 



four-quark operator, just as in the case of B — B mixing. The [AS"! = 1 and [AS*! = 2 Wilson 
coefficients of this effective field theory are evolved down to the scale /i^c = 0{mc) at which the 
second step of the calculation is performed: Now the bottom and charm quarks are integrated out and 
the effective field theory set up in the first step is matched to another effective field theory. The new 
theory treats and nic as heavy scales, so that all box diagrams involving at least one charm quark 
are effectively contracted to a point. All information on nic (and mb which plays a minor role) resides 
in the Wilson coefficient of the local AS" = 2 operator 

Q = di'luSLdiYsL- (121) 

The effective |A5| = 2 Hamiltonian can therefore be written in a similar way as the |Ai?| = 2 
Hamiltonian of Eq. dlOll) : 

^1A5|=2 ^ [{VtsVCafvttSixt) + 2VtsV:^VMv,tS{x,,xt) 

+ {VcsV*^fr]ccXc] hK{pK)Q{l^K) + h.c. (122) 
The NLL results for the short-distance QCD factors read 

(1 3 Cjc \ 
• (123) 
rric J 

The QCD coefficients in Eq. (11231) were calculated to LL accuracy in Ref. [30]. The NLL calculation 
of r]tt [18] is analogous to that of rjB, with one new feature: When crossing the threshold fi^c one 
must change the number of active flavours in the QCD (3 function and the NLL anomalous dimension 
7+ from / = 5 to / = 3. The NLL results for rjd [31] and 7]cc [32] have a sizable uncertainty, 
because they are sensitive to the low scale of ~ rric where ag is large, rjcc also exhibits a sizable 
dependence on as{Mz) and on rric = mc{mc), so that the central values quoted in the literature vary 
over some range. The expression in Eq. (11231 ) approximates the dependence on ttIc and corresponds 
to as{Mz) = 0.119 lb 0.002. The scale /xa' must be chosen below rric and is typically taken around 
IGeV, where perturbation theory is still applicable. One finds bxilJ'K = IGeV) = 1.24 it 0.02 and 
the error stems from the uncertainty in a^. 
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In the discussion of |A5| = 2 transitions we must also address corrections of order "i-^j^j^j/m^^j^^y 
which correspond to subleading terms in the operator product expansion of Eq. dTQ] ). While these 
corrections are of order Aq^j-j/m^ for the first term in ijl^'^N^^ jj^gy ^j-g gf order Aq^^j-j/m^ in the 
case of the charm contributions involving S{xc,xt) = 0{xclnxc) and S{xc) ^ Xc in Eq. (I122I ). 
The largest of these power corrections involves two |A5| = 1 operators and corresponds to the box 
diagram in Fig. [T] with two internal up-quarks. To understand the power counting, recall that the charm 
contribution in //l^-^'N^ is proportional to M^^Xc = m^, while the box with up-quarks involves no 
power of TTic, so that its size is characterised by the hadronic energy scale Aqcd- Including this bilocal 
contribution we write: 

Mi2 = ^— {K\h\^^\=^\K) - Disp -!— [ d^x {K\h\^^^=\x) h\^^\=\0)\K) . (124) 
2mK ^rriK J 

Here "Disp" denotes the dispersive part of the matrix element, which is introduced in Eq. (I24l) and is 
discussed after Eq. (11041 ). The enhancement of the second term stems from the so-called A/ = rule 
which describes the non-perturbative enhancement of the decay Kshort (vr7r)/=o- The two terms in 
Eq. (11241 ) are usually referred to as short-distance and long-distance contributions. The long-distance 
contribution has defied any reliable calculation from first principles so far. In this humbling situation 
we can only compare the experimental value of AMk to the short-distance contribution 

AMlP = K^^^^t^. (125) 

rriK 

In order to compute AM|P we need the hadronic matrix element 

{K\QM\K) = \Mlfl-^ (126) 

Contrary to the situation in the B system, the Kaon decay constant fx = 160 MeV is well-measured. 
We remark here that we know Bk in a particular limit of QCD: If the number of colours Nc is taken to 
infinity, {K\Q{iik)\K) can be expressed in terms of the current matrix element {Q\dL^i,SL\K) which 
defines Jk- For Nc = oo one finds Bk /bxilJ-K) = 3/4; including certain calculable ("factorisable") 
l/Nc corrections changes this to Bk /bxifJ-K) = 1. A recent lattice calculation finds [33] 

Bk = 0.72 ±0.04. (127) 

The experimental value of the i^^iong-^short mass difference is [8] 

AM^'^P = (3.483 ± 0.006) ^eV = (5.292 ±0.009) • 10"^ ps"\ (128) 

Inserting Eqs. (11221 ) and (11261 ) into Eq. (11251 ) gives 

AM^^ 

= (0.98±0.22)i?^. (129) 

K 

AM^ is dominated by the term proportional to (VcsV*^)"^ and the error in Eq. ( 11291 ) essentially stems 
from rjcc in Eq. (11231 ). This uncertainty will shrink when rjcc is calculated to NNLL accuracy. With 
Eq. (11271 ) we find that ijl^'5'l=2 contributes (70 ± 25)% to the measured AMk- 
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The off-diagonal element of the decay matrix is given by 

= Abs^— [ d'^x {K\H\^^'^\=\x) H\^^\=\0)\K) (130) 
2mx J 

= ^E(2^)''^'(^'^-^^/)(^l^'^'"=V>(/l^l^''l='l^) ^ ^^oAo. (131) 

Here "Abs" denotes the absorptive part of the matrix element. is an inclusive quantity built out of 
all final states / into which both K and K can decay. A special feature of the neutral Kaon system is 
the saturation of by a single decay mode, which is K {tt-k)j=q. The notation Aq and Aq for the 
corresponding decay amplitudes has been introduced after Eq. (fT4l) . is a non-perturbative quantity 
and its computation on the lattice involves the difficult task to understand and master the AI = rule. 
The relation between and AFk has been derived in Eq. (l45l) . Experimentally we have [8] 

Ar^P = (7.335 ±0.004) /icV = (11.144 ±0.006) • 10"^ ps"^ (132) 

With Eqs. (11281 ) and (11321 ) we have precise experimental information on |Afi2| ~ AMk/2 and 
|ri2| ^ Arx/2. To fully characterise K — K mixing we also need to know the phase cj) defined 
in Eq. (|38] ). As in the case of B — B mixing we study a CP asymmetry in a flavour-specific decay 
mode. With Eqs. (ITTI ) and (|4T]) one easily finds 



Al = 



r(Kiong ^ i+l^TT-) + r(i^iong ^ ^-l>7^^ 

l + \q/p\^ ^ 2- 



(133) 



At this point it is worthwhile to look back at the quantity ex which we have encountered in the first 
lecture in Eq. (fT4l) . From Eq. ([TTl) we have learned that Re ex measures CP violation in mixing 
quantified by 1 — \q/p\, just as Al in Eq. (11331 ). While Imex is related to a different physical 
phenomenon, namely mixing-induced CP violation, it provides the very same information on the 
fundamental parameters of K — K mixing: Since K (7r7r)/=o dominates ri2, the CP-violating 
phase of Aq/Aq equals arg ri2, see Eq. (11311 ). With this observation and the help of Eq. (|42)) we can 
express ex in Eq. ([17] ) entirely in terms of AMk, ATk and (p. Interestingly, the phase of ex (see 
Eq. ([14]) ) is simply given by 

= arcan^, (134) 

Al K/2 

More details of this calculation can be found in Chapter 1.6 of Ref. [34]. Nature chose AMk ~ 
Arx/2 by accident, so that (p^ in Eq. ([14] ) is close to 45°. The bottom line is that (p^ carries no 
information on CP violation and that | \ and A^ involve the same fundamental CP- violating quan- 
tity, which is 0. To extract (p from ^4^ in Eq. (1133b or from ex in Eq. ([TT] ) we use Eq. ([42] ). with 

2|Mi2|/|ri2| ~ AMK/{ArK/2) traded for tan 

Al = ^sin(2</),)0 + O(</.2) 

eK ^ ^sm{(l,,)e'^^<P + 0{(l)^) (135) 



Using the experimental value 



A^^P = (3.32 ± 0.06) X 10" 
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gives 



(6.77 ±0.12) X 10" 



(136) 



This number is in reasonable agreement with = (6.48 ±0.03) x 10 ^ found from ek with Eq. (I135I ). 
Next we relate (pto a. constraint on {p,ri): Specifying to the standard phase convention for the CKM 



matrix (with VusV*, real and positive) we start from Eq. (1381 ) to write 



arg 



M 



12 



12 



where 



ImMi2 

|Mi2| 



arg(-ri2) 



arg(— r 



12 



ImMi2 



±e 



K 



arg 



(137) 



(138) 



In Eq. (11371 ) I have used that the phases of M12 and — ri2 are separately small in the adopted phase 
convention and further traded |Afi2| for the experimental AMk/2. In Eq. (11381 ) the saturation of ri2 
by ^0 Aq in Eq. (11311 ) has been used. Thus —^k is just the CP-odd phase in the decay K (7r7r)/=o- 
A recent analysis has estimated ~ —1.7 • 10~^ [35], so that contributes roughly —6% to the 
measured value of 4>. The dominant term proportional to Im M12 = Im {K\H^'^^^^'^\K) involves the 
CKM factors 



-lm{VMf: 



P) 



(139) 



where the lowest-order Wolfenstein expansion has been used. Inspecting the dependences of the CKM 
factors on p and rj we see that the experimental constraint from (j) defines a hyperbola in the (p, fj) 
plane. Combining Eq. (11391 ) with Eqs. (|122l) and (I137I ). inserting the QCD factors from Eq. (11231) and 
the matrix element from Eq. (11261) and finally using (f) = (6.48 ±0.03) X 10-3 from ex this hyperbola 
reads 



V 



1 0.34 ± 0.03 
1^ 1.3±0.1-p' 



(140) 



The uncertainties in Bk from Eq. (11271) and from rjcc and rjct in Eq. (11231) (reflected by 1.3 ± 0.1) 
inflict errors of similar size on the rj extracted from Eq. (11401 ). The numerator 0.34 ± 0.03 is calculated 
with \ Vcb\ = AX'^ = 0.0412 ± 0.0011. The 10% uncertainty of this number stems solely from the 
error in \ Vcb\, which enters rj in Eq. (11401) with the fourth power. 

The neutral Kaon system is the only neutral meson system for which all three quantities AM, AT 
and (j) are measured. It should be stressed that also the sign of AT/ AM is firmly established. Mea- 
suring sign (AT/ AM) is difficult for all meson-antimeson systems. In the neutral Kaon system the 
measurement of AM and sign (AF/AAf) uses Kshort regeneration: If a Kiong beam hits a nucleus in 
a target (the regenerator), strong inelastic scattering changes the |i^iong) state into a superposition of 
l^iong) and |Kshort) giving access to observables which are sensitive to AM and the abovementioned 
sign. For details on these experimental aspects I refer to [36]. 

Finally I discuss D — D mixing: Box diagrams in Fig. [T] with one or two internal b quarks are highly 
CKM-suppressed. The dominant box diagrams with internal d and s quarks suffer from a very efficient 
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Figure 7: Gluonic penguin diagram with an internal top quark. 



GIM suppression proportional to m^/m^. This makes the diagrams sensitive to very low scales and 
perturbative calculations of AMd, ATd and are put into doubt. In the effective theory both M12 
and are dominated by the bilocal contribution with ijl^'^l=^. The only possible clear prediction 
is the qualitative statement that all these quantities are very small. Theoretical calculations usually 
quote numbers for the quantities x = AMd/Td and y = Ar£)/{2T£)). The theoretical predictions 
for l^l, \y\ cover the range from zero to \y\ ~ 0.01 and come without reliable error estimates. 
Therefore current experimental values are compatible with the SM but may also be dominated by a 
new physics contribution. A "smoking gun" of new physics, however, would be the discovery of a 
non-zero CP asymmetry in the D system. 

3.4 Mixing-induced CP asymmetries 

At the end of Sect. 12.31 we have learned that mixing- 
induced CP asymmetries can provide clean information 
on fundamental CP phases in the Lagrangian. These 
CP asymmetries involve the interference between mix- 
ing and decay amplitudes as depicted on the right. 

In this lecture we restrict the discussion to gold-plated modes which involve a CP eigenstate fcp in the 
final state, cf. Eq. (fTOtFl In the Bd and Bg meson systems the mixing-induced CP asymmetries are a 
real gold mine, because there are many decay modes satisfying the condition for a golden decay mode 
as defined after Eq. (l62l ). Prominent examples are the decays Bs —>■ J/ipcl) and B^ —>■ J/ipKshovt, 
whose decay amplitudes essentially only involve the CKM factor VcsV*^^. To understand this first 
note that the decay proceeds at tree-level by exchanging a W boson. There are also contributions 
involving an up, charm or top quark loop, with attached gluons splitting into the charm-anticharm pair 
hadronising into the J/ip meson. Such diagrams are called penguin diagrams. A penguin diagram in 
the narrow sense only involves one neutral vector boson (which can be a gluon, photon or Z boson). 
A gluonic penguin diagram is depicted in Fig.|7] (Yet a. J/ip cannot be produced from a single gluon. 
One needs a photon or three gluons at least.) In the context of mixing-induced CP asymmetries 
one often speaks of penguin pollution, because the penguin diagrams may involve different CKM 
factors than the tree diagram spoiling the golden-mode property. To estimate the penguin pollution in 
Bs — > J/i/jcl) and Bd J/tpKshovt first use the unitarity relation VtsV^l = —VcsV*f^ — VusV*^^ to write 

^|AB|=i ^ VcsV:,h, + V^sVcbhl + VusV:,hu + CKfe/il- (141) 

Here the last two terms are highly suppressed, since ~ 0.03 Moreover, hu has no 

tree contributions, but solely stems from penguin diagrams with up and top quarks. Since these loop 

One can also identify gold-plated decays into non-CP eigenstates, important channels are e.g. 



B 



q/p 



B 



At 
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effects involve non-perturbative physics, it is difficult to quantify the loop suppression. Still, the CKM 
suppression is efficient enough to render the modes gold-plated at the level of a few percent. Since the 
CKM elements are factored out in Eq. (I141I ). he and hu only contain Wilson coefficients, operators 
and real constants. Importantly, /i^ c and /i^ c are related by the CP transformation: 

hi,c = {CP)hu,cCP (142) 

While I discuss Bs J/'4'4' and B^^ J/ipKshort here for definiteness, the results apply to other 
gold-plated M — > /cp modes as well, with obvious replacements for the CKM elements. The under- 
lying reason for the cancellation of hadronic uncertainties in gold-plated decays is the CP invariance 
of QCD: While we cannot compute (/cp|^c|^fl we can relate this matrix element to {fcp\hl\B) 
through 

{fcp\hijB) = {fcp\{CP)^hu,cCP\B) = -vcp if Cp\hu,c\B), (143) 

where I just used the CP transformations of Eqs. (l9l-IT0l) and Eq. (11421) . We first apply this to the 
decay mode Bs J/tpcl). The final state consists of two vector mesons. By conservation of angular 
momentum they can be in states with orbital angular momentum quantum numbers / = 0, 1 or 2: The 
two spin-1 states of the vector mesons can be added to a state of total spin 0, 1 or 2, which requires an 
orbital angular momentum of Z = 0, 1 or 2 to give a J/ipcj) state with zero total angular momentum. 
The p-wave state with / = 1 is CP-odd and the other two states are CP-even, owing to the parity 
quantum number (—1)' of their spatial wave function. Experimentally one separates these states by 
an angular analysis [37, 38] of the data sample. This can be done including the full time dependence 
of the decay, so that we can isolate the time-dependent CP asymmetries in the different partial-wave 
channels. The most-populated state is the CP-even I = (i.e. s-wave) state. Writing fcp = {J/t(jcl))i 
with 7]cp = {—ly we obtain for the amplitudes Aj^p and Af^p (see Eq. (fT2])): 

^/CP ^ (/cp|gl^^l=^|^s) ^ V:sVcb{fcp\hl\Bs) ^ 

Af^p (^f^p\H\m=^Bs) v,sV:,{fcp\K\Bs) '"'Ve.y;, 

Combining this result with Eqs. (|47] ) and (fT3] ) we find 

A/CP = VCP^^ = .cPe-/^^ (145) 

In the last step I have used the definition of (3s in Eq. (177] ). With Eq. (11451 ) we can calculate the 
time-dependent CP asymmetry of Eq. ( [6T| ). First we verify that our golden mode satisfies | A j^p | = 1, 
so that A'^p in Eq. (|62l ) vanishes. The other two quantities in Eq. (|6T| ) evaluate with Eq. ( 11451 ) to 

^CP ~ ~^CP sin(2/3s) and A^r = —Vcp cos(2/3s), so that (neglecting the tiny 0{a) term) 

sm(2l3s)sm(AMB t) 
cosh(ArB, t/2) - ricp cos{2f3s) smh^ATs^ t/2) 

In the SM f5s is small and ai^jj^^-^^ (t) is an ideal testing ground to find new physics [27, 38]. 

Next I discuss Bd J/t/jK^hort- The final state has orbital angular momentum / = 1 balancing 
the spin of the J/^p. Neglecting the small CP violation m K — K mixing we can regard the i^short 



*In flavour physics matrix elements like ^ | A^) are always understood to include the strong interaction. This 

means that the fields are understood as interacting fields in the Heisenberg picture with respect to the strong interaction. 
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Figure 8: Interfering amplitudes which give rise to mixing-induced CP violation for the two golden 
modes discussed in the text. 



as CP-even. The J/t/j is CP-even as well and the orbital angular momentum contributes a factor 
of —1 to the total CP quantum number. Thus Vj/ipK^hort ~ From Fig. [8] we observe a novel 
feature compared to Bg —^ J/ifjcj). The interference of the and Bd decays involves K—K mixing: 
The Bii decay involves the K component of i^short. while the Bd decays into the K component of 
ETshort- Experimentally the Kshort is detected via a pair of charged pions whose invariant mass equals 
Mk, denoted here by {tt^'k^)k- Therefore we should identify the amplitudes ^/cp=J/V'i^short ^'^^ 
^/cp=j/,/.i^,hon with A{Bd JI^K J/V'(7r+7r-)/^) and :4(5d ^ Jj^K J/V'(vr+7r-)i^), 
respectively. Therefore 

In the last step I have used the definition of /? in Eq. (1721 ) and neglected arg[— 14^1/* ~ 
A^X'^l) < 10~^, so that Im Aj/^^-^^^^^ ~ sin(2/3). We may further neglect AT^^ in Eq. (|6TI ) to find 
the most famous time-dependent CP asymmetry, 

« J/Vi^shon (t) = sin(2/?) sin( AMe.t) . (148) 



Finally I give a (very incomplete) list of other golden M —>■ fcp decays. The decay Bs —>■ J/i/^cl) 
can be substituted for Bg —> J/ipr]^'\ which does not require any angular decomposition. In a hadron 
collider experiment r/'s and ri"s are hard to detect, but Bg J/'ipij^'^ is interesting for B factories 
running on the T{5S) resonance. While the modes discussed above provide insight into the physics of 
B—B mixing, one can also use mixing-induced CP violation to probe CP phases from new physics in 
loop-induced B decays such as Bd —>■ 0Kshort [39]. This mode is triggered by the quark decay b — > 
sss. The same transition in probed in Bs — > 4'4'. Likewise new physics in the b — > sdd amplitude may 
reveal itself in Bs — > i^short^short- Gold-plated decays are KshortT^^ and KshovtP^, 
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wiiich are penguin-free c — > sdu decays. A gold-plated K decay is Kiong tt^vV [40]. Here no 
meson-antimeson oscillations are present, but K — K mixing nevertheless enters the process through 
the mass eigenstate A'long- The final state tt^vV is CP-even and the dominant contribution to the 
decay involves mixing-induced CP violation, i.e. the decay amplitude is proportional to ImAj (see 
e.g. Ref [41]). 



3.5 The unitarity triangle 

Many measurements contribute to the global fit of the unitarity triangle defined in Eq. (TTOI ) and de- 
picted in Fig. [3] Conceptually it is useful to disentangle tree decays from FCNC processes: Tree-level 
ampUtudes are insensitive to new physics and therefore determine the true apex (p, 77) of the unitarity 
triangle. In principle one could determine the unitarity triangle in this way, insert the result into the 
SM predictions of the FCNC processes and then assess the possible impact of new physics on the 
latter. In practice, however, the tree constraints still suffer from large uncertainties, while for example 
aj/^ft:,hort(^) ^1- <I148I ) and AMb^^/AMb, in Eq. (11141) determine sin(2/3) and the side Rt (see 
Eq. (ItT]) ) fairly precisely. Therefore, for the time being, it is best to combine all information into a 
global fit of the unitarity triangle. 

From b — > ciu decays \ Vcb\ — AX^ is precisely determined. Therefore we realise from Eqs. (T/T]) and 
(1761) that any measurement of \Vub\ essentially fixes the side i?„ of the triangle. \Vub\ is determined 
from (inclusive or exclusive) semileptonic b ^ u decays and hadronic uncertainties limit the accuracy 
of the extracted \Vub\ to 8-10%. The theoretical methods used to determine \Vcb\ and \Vub\ are briefly 
reviewed in Ref. [42]. The angle 7 of the unitarity triangle is currently measured in two ways from 
tree-level decays: First, the interference of the b cus and b — > ucs amplitudes in — > ''D ^K^ 
decays is exploited [43]. Second, one measures mixing-induced CP violation in Bd — > vrvr, Bd pir 
or Bd pp decays, which allows to find the angle a of the desired triangle. These modes are 
not gold-plated and suffer from penguin pollution, which, however, can be eliminated by means of 
an isospin analysis [44]. While the extracted result for a is sensitive to new physics in Bd — Bd 
mixing, this possible effect can be eliminated if the measured a™P and /3™p are combined to give 
^exp = _ Q,cxp _ ^cxp Combining the constraints from \ Vub\, 7 and a with those from meson- 
antimeson mixing discussed in this lecture results in the unitarity triangle shown in Fig.|9] 



Suggestions for further reading 

There are many good review articles on meson-antimeson mixing and flavour physics in general, 
putting emphasis on different aspects of the field. A student interested in the theoretical foundation 
of flavour physics, effective Hamiltonians and higher-order calculations is referred to the lecture in 
Ref. [48] and the review articles in Refs. [34,49]. Most reviews and lectures focus on CP violation 
and I recommend Refs. [47] and [50]. I have only briefly touched D — D mixing, two review articles 
dedicated to D physics are cited in Ref. [51]. Lectures covering both K and D physics can be found 
in Ref. [52]. A concise summary of the physics entering CKM metrology can be found in Ref. [42], a 
more elaborate article on the subject is Ref. [53]. Standard textbooks on flavour physics are listed in 
Ref. [54]. 
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